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Abstract 

We presented a general classification scheme for the tripartite L x N x N en- 
tangled system under stochastic local operation and classical communication. The 
whole classification procedure consists of two correlated parts: the simultaneous 
similarity transformation of a commuting matrix pair into a canonical form and 
the study of internal symmetry of parameters in the canonical form. Based on this 
scheme, a concrete example of entanglement classification for a 3 x x system 
is given. 

PACS numbers: 03.67.Mn, 03.65.Ud, 02.10.Xm 



1 Introduction 

Entanglement represents the essence of quantum theory compared to the classical 
theories. With the development of quantum information science, entanglement is now 
expected to be generated from a wide range of systems and is generally regarded as 
the key physical resource in implementing various quantum information tasks pQ. The 
classification of entanglement is generally a mathematically difficult problem. It has been 
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pointed out that only in a kind of tripartite system, which has a qubit as one of the 
subsystems (2 x M x A^), one may find finite number of entanglement classes [2]. The 
complete entanglement classification of a 2 x M x pure system under stochastic local 
operation and classical communication (SLOCC) has been done in Refs. [3], H], where 
canonical forms are constructed via Jordan decomposition of the matrices. There the 
eigenvalues of the Jordan forms serve as the continuous parameters in the entanglement 
classes [3]. Despite the great progress made in 2 x M x systems, in many classification 
schemes the classification of 3 x x turns out to be a "wild" problem [5]. Attempts 
toward the classification of a multipartite system in literature mainly concentrate on its 
tensor ranks or local ranks of the tensor form of the quantum state, and there is still 
no systematic method for constructing the canonical form of L x N x N under SLOCC 
in literature to our knowledge. Recently, a method featured by the decomposition of 
the entanglement classification under local unitary transformations into different parts 
has been proposed [6J. In this work, by decomposing entanglement classification into 
two complementary ingredients, we propose an entanglement classification scheme for a 
tripartite L x N x N system under SLOCC. 

The structure of the paper is as follows. In Sec. 2, by representing the L x N x N 
quantum state in the form of a three-way tensor, its entanglement classification under 
SLOCC is transformed into two correlated problems: the simultaneous similarity trans- 
formation of a commutating matrix pair into a canonical form and the representation of 
the symmetry of parameters in the canonical form. In Sec. 3, a concrete example of 
entanglement classification for a 3 x A^ x A^ system is given in which the explicit form of 
the symmetry property of the state is shown. A brief summary and concluding remarks 
are given in Sec. 4. 
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2 The classification scheme 



Considering a general tripartite pure quantum state with dimensions L x N x N, the 
wave function can be written as 

L,N,N 

^= E (1) 

i,j,k=l 

where \j), \ k) are bases of the three parts respectively. The complex coefficients 'jijk 
can be treated as a three-order tensor which then turns into the form of a tuple of L 
complex N X N matrices, i.e., 

^ = (ri,r2,---,ro, (2) 

where Fj, i G {1, 2, ■ ■ ■ , L} are N x N complex matrices and have [jiljk as their elements. 
A state ip' = (r'^, r'g, ■ ■ ■ , r^) is said to be SLOCC equivalent [2] to '0 in the case of 

L 

J^T.^PT, Q = T[,t = {1,2,. --,1}. (3) 

i=i 

Here, T (where its elements are Tjj) is an L x L invertible matrix, i.e., invertible local 
operator (ILO) acting on the first partite, and P and Q are N x N ILOs acting on the 
second and third partite, respectively. 

For a tuple of matrices (Fi, F2, ■ ■ ■ , F^), there always exists a tuple of numbers {ti, t2, - ■ ■ , ti) 
ti G C, such that the linear combination of FjS, 

L 

F = 5^t,F, (4) 

gives the maximum rank, the r(F). (We denote the rank of a matrix as r(-) hereafter.) 
For every quantum state in the form of Eq.([2]), we can always transform them into the 
form where the Fi has the maximum rank by performing the transformation T, of which 
the first row is chosen to be the tuple of (^1,^2, ■ ■ ■ ^^l), i-e., Tij = tj. Therefore, in the 
following we can then restrict our discussion of the quantum state classification to that 
of the tuple (Fi, F2, ■ ■ ■ , F^), where Fi has possessed the maximum rank. 
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The strategy we adopt in the entanglement classification in this work is quite trans- 
parent, which can be roughly outlined as the manipulation of three commutative ILOs 
T, P, Q in different order and combination. We first transform the quantum state in a 
tuple of matrices according to specific rank property for each matrix by a subset of ILO 
T, then obtain a canonical form by using ILOs P, Q, and finally achieve the entangle- 
ment classification by reconsidering the rest a subset of ILO T as a degeneracy condition 
between the canonical forms. Following, we explain in detail the classification procedures, 
in which two situations of matrix Fi with full and nonfull ranks are discussed separately. 

2.1 The case of Ti with full rank 

In this case, obviously every quantum state in form ([2]) can be transformed into the 
following form by ILOs P, Q defined in Eq.Q 

ij = {E,T2,--- ^Tl) , (5) 

where i? is a unit matrix. Then a subsequent ILO transformation 

/ 1 ■■■ o\ 

-t21 1 ■■■ 

T= : : ... : ' (6) 

\-tLl ■■■ 1/ 

where tn is one of the eigenvalues of Fj, can always be applied to -0 and makes r(Fj) < 
after the transformation. By definition, two such quantum states are equivalent 

under the transformations of ILOs P, Q if and only if 

{PEQ, PT2Q, ■ ■ ■ , PTlQ) = {E, F'2, ■■■ Xl), (7) 

which tells Q = P^^. Equation ([7]) now turns into a simultaneous similarity relation of 
matrices (F2, ■ ■ ■ , F^), 

{PEP-\ PT,p-\ ■■■ , PTlP-') = {E, F'2, ■ ■ ■ , F'J . (8) 

According to ff, '8] , all these kinds of problems can be transformed into the simultaneous 
similarity of a pair of commuting matrices, i.e., Eq.® is equivalent to 

{E, PA2P-\ PA,P-') = {E, A'„ A',) , (9) 



where ^2,3 are composed by Fj, and [742,^3] = 0. This makes the entanglement classifi- 
cation of a 3 X X system of great importance. The canonical form of Eq.Q under 
the similarity transformation can be written as {E, J, A}, where J is the Jordan form of 
A2, A is the canonical form of ^3, and [A, J] = 0. J can be represented as J = ^ JnX\)i 
where Jni{\) are rij x rij matrices which take the following form: 








1 

A. 





1 

A,; 







0\ 






1 

A,; 



(10) 



Due to [v4, J] = 0, the matrix A generally has the form oi A = @ Ai accordingly, where 
Ai is a nj X Ui triangular Toeplitz matrix (see Appendix Rl) : 



A,- 





an 


ai2 


ais ■ 









an 


ai2 ■ 


aini-2 








aio 


an ■ 


aim— 3 











aio ■ 


aim— 4 










■ 


aio J 



111 



Since the Jordan block can be defined as a polynomial function of J„i(0), 

Jm{^^) = Ai + J„,(0) , 
the matrix Ai can also be similarly defined as 

rii-l 



;i2) 



A=5^a.„<(0) 



(13) 



n=0 



By defining f[A,]{x) = En=o ^i"^"' ^e have Ai = f[A,]{Jmi^))- 

However, there still exist other transformation Ts that can relate i/j = {E, r2, ■ ■ ■ , Tl) 
to another quantum state ip' = {E,r2,- ■■ ,r^). The effect of the transformation T is 
introducing a degeneracy on the canonical form of {E, J, A). That is, the states {E, J, A) 
and {E, J', A') with parametric difference may belong to the same entanglement class. 
Define a set a = {Aj, a^lO < n < rii — 1}, which contains all the parameters in {E, J, A). 
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The degeneracy can then be expressed as a symmetry transformation V on the parameters 
of the canonical form, i.e., 



a' = Pa . 



(14) 



Thus, the canonical form is of the entanglement classification of the quantum state up to 
the symmetry V. 

2.2 The case of Fi with non-full rank 



In this case, all of the quantum states can be transformed into the following form by 
ILOs F, Q 



^ = (A,r2,--- ji 



(15) 



Here, A has the following form: 



A 



E(N-i)x{N-i) 







0,; 



/I •• 


• 


■■ 


■ 0\ 


■■ 


■ 1 


■■ 


■ 


■■ 


■ 


■■ 


■ 


Vo ■■ 


■ 


■■ 


■ 0/ 



(16) 



with i > . 

Then we apply another type of ILOs, Pa,Qa, such that they keep A invariant and 
transform all the Tj, j G {2, ■ ■ ■ , L} to a form of direct sum of two matrices. The 
philosophy here is the same as that of the less than full rank case in [3], that is. 



/Till ■ 


■ ■ 7jln 


■ 


■ \ 


Tjnl ■ 




■ 





■ 





/?.ll ■ 


Pjlm 


V ■ 










(17) 



Here, N = n + m, and the matrices {3j cannot be simultaneously for j and cannot be 
further decomposed into the direct sum of submatrices via Pa,Qa- From this and from 
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the condition that the first matrix A has the maximum rank in comparison to Tj, we can 
also deduce m > i. With the partition (|T7|) . A can be partitioned accordingly as 



A 



En 



A' 



/I •• 


■ 


■■ 


■ 0\ 


■■ 


■ 1 


■■ 


■ 


■■ 


■ 


1 ■■ 


■ 


Vo •• 


• 


•• 


• 0/ 



where 



A' 



E (m—i) X (m—i) 







0,; 



(19) 



Now from the direct sum forms of Eqs. (fT7l [T8l) . the task of classification under Pa, Qa 
can be divided into two problems: the simultaneous similarity transformation of matrices 
7j and the construction of matrices Pj in canonical form. The first one can be solved 
with the help of results in Sec. 12.11 For the second problem the necessary and sufficient 
condition for the canonical form of (3j goes as follows: 

L 

rmax(A' + ^ ajl3j) = m - i , aj e C . (20) 

Similar to the full rank situation, the transformation T here will also induce a sym- 
metry V to the canonical form, which can be schematically expressed as 



(21) 



where b = (Aj, a^^, /3jh)- 



In all one can conclude that for each matrix triple {E,J,A), the following theorem 
holds: 



Theorem 2.1 The classification of Lx N x N can be represented in a triple-matrix form 
{E, J, A) up to a parametric symmetry V. 
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Now we note that the entanglement classification of quantum state L x N x N under 
SLOCC can be decomposed into two problems: the simultaneous similarity transformation 
of a commutating matrix pair and the symmetry analysis of the corresponding canoni- 
cal form. Therefore the difficulty in entanglement classification then turns to that in the 
symmetry property study, which requires further analysis for different entanglement struc- 
tures. Analysis of the symmetry properties represents the essence of the entanglement 
class in our classification scheme. In the following we show the explicit form of symmetry 
in the classification of the 3 x N x N quantum system, from which the symmetry nature 
of state L X N X N can be figured out. 

3 Entanglement classification of 3 x x system 

In this section we give an explicit example of a 3 x x system: i/j = [Ai, A2, A3), 
where r{Ai) = N (full-rank case), and when transformed into {E, A2, A'^) via P, Q it has 
[A2, A'^] = 0. From the discussion in Sec. 2 it is clear that the canonical form of this state 
is {E,J,A), where J is the Jordan form of A2 and [J, A] = 0. We can always achieve 
r{E) > r{J) > r{A) by applying a certain transformation T. In order to realize the 
unique entanglement classification, we need to find out the specific form of symmetry V. 

An invertible matrix T can be decomposed as T = FLDU, where L {U) is the lower 
(upper) triangular with all diagonal entries equal to 1, i5 is a nonsingular diagonal matrix, 
and P is a permutation matrix [9]. Among these four matrices F, L, D,U: (1) P and L 
are now contradicting with r{E) > r{J) > r{A) in most cases, except for very limited 
(less than the number of distinct values of Aj) and isolated values which can be treated 
individually; and (2) D corresponds to a trivial rescale of the whole matrix which can 
be treated by the renormalization of quantum state. Now the remaining effective T 
transformation is 

/I ^1 /I ^1 0\ /I /I 0\ 

[7= 1 Z3 = 1 1 1 ^3 - (22) 
\0 1 / \0 1/ \0 1 / \0 1 / 

The admissible T matrix which can be applied on the canonical form {E, J, A) is now 
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decomposed into a sequence of three operations, like 
'E + ziJ\ fE\ ^(.2) fE + Z2A 




J 
A 




(^2) 



J 

A 




\.JA] 



E 



> I J + ^3A I , (23) 



where T^-^y] is an elementary operation that multiplies Y with z and adds it to the X 
column. Because all the E, J, A are direct sums of block-diagonalized submatrices and 
have the same partitions, we can therefore treat Eg. (1231) with their diagonalized blocks. 



3.1 The first superposition T, 



(^i) 

[EJ] 



After the transformation Tj^Jj] , the initial canonical form turns to 




(24) 



which is no longer a canonical form. To obtain the canonical form corresponding to it, 
operators P and Q should satisfy 



PiE + zM\))Q = E 
for each block. This equation can be reexpressed as 



(25) 



(26) 



where 5'[Bj](a;) = 1 + ziXi + zix. The reciprocal of polynomial function g[^]}{x) can be 
readily obtained, that is (note J^-(O) = if n > rij). 



Cfel(^n.(0))= 9\l\\{JnM) 



rii 



1 + ZiXi ^ VI + ^l^i 
n=0 



(27) 



Since Eq.(l26D leads to F = Q-^ ■ Cfe{(Jn,(0)), we can further get 

PJnA^i)Q = Q~'-Cfei(-^n.(0))-^(A.)-Q 



A,; 



1 + ziXi 
•(-i) , 

[(BJ)J]> 



^§(l + ^iA.)^(l 



+ ziXi 



n-l 



4"(o) 



■Q 



Q'' ■ fliSjYjpnM) ■ Q ■ 



{2i 



Here /[^(Ej)j](^n»(0)) is defined to be C[ej] (-'nilO)) ■ J„,-(Aj), which represents the influence 
of the superposition of E and J on J after the operation P, Q. Similarly, we can obtain 
this operation influence on Ai, i.e., 



PAQ 



P ■ flA.]{Jn..m ■ Q 



_n=0 
-1 ^{^i) 



^ 1 + ziXi 5^ ( 1 

n=0 m=0 ^ 



-Zl 



Q 



(29) 



In order to see what the canonical form becomes now, we perform a similarity transfor- 
mation S on the second matrix, Eq. (l28|) . and transform it to the Jordan form 

A, 



which can be verified with the help of ( l28i) . From ( l30l) we then have 



(30) 



S-'Q-'- J„,(0)- QS = f, 



KEJ)J] 



J, 



rii 



1 + ziXi 



(31) 



Here f[^Ej)j] inverse function of f^^Ej)j]- This similarity transformation will also 

make 



(^i) 



o J„ 



[(EJ)A] ^^[(i?J)J] '^"^^\l_^z^X 



A,; 



(32) 



where the Jordan canonical form is treated as a polynomial function [see Eq. (fT2l) ] and 
f°9{x) = figix)). 



We note that when the T transformation of E + ziJ applies on 

iE,JnAk),A,) , 

where Ai = /[yi.](J„. (0)), it will make the canonical form transform into 

A,; 



(33) 



{E, Jn 



1 + ^lA ' ' ' 
10 



(34) 



where A'- = f[^Ej)A] ° f[{Ej)j] ° -^n, (yT^^^-} ' '^'^^^^h dM]) and are similar in 
block partitions with different parameters, they are equivalent entanglement classes. That 
means the entanglement classes parameterized by a = (Aj, aj„) are not inequivalent with 
different parametric values, while parameters undergo a transformation Vz^ characterized 
by zi, i.e., 

A,, 



A,; 



1 + ^1 A. ^ ^35) 



a'- 



where is the entries of A'- defined in Eq. (l511) . or simply, 

(36) 















A') 




u 



3.2 The second superposition Tj^J^j 

Similar to the above section, the transformation requirement 

P{E + Z2Ai)Q = E (37) 

leads to 

P- glEilJuM]- Q = E, (38) 

where the polynomial glEA]{x) = X]n=o^"^"' "with uq = 1 + Z2aio, aj = Z2aij when j > 0. 
Define C!^i{JnXO)) = [^?Sai(^(0))]"^ such that P = ■ Cl^iiJnX^)). What we need 
to know now are the values of PJn^{\i)Q and PAiQ. 

For P J„.(Ai) (5, we have 

PJn,(A.)Q = Q-^ ■ Cfci(Jn,(0)) ■ ^(A,) ■ Q 

= Q'' ■ f[iE\)J]iJn.m ■ Q . (39) 



11 



Here /gi](Jn,(0)) = &Jn,m ■ JnX\)- For PA,Q, we have 

PA,Q = P- f[A,]{JnM)-Q 

There always exists a similarity transformation such that 

S-'Q-'ft^%j^{JnM)QS = Jn, ^^^^ 



and hence 



Then we have 



^ 1 + Z2aio 



(40) 



(42) 



s-'Q-'flli\^4JnM)QS 

= fliEA)A] (S-'Q-'JnMQS) 

We note that when the T transformation of + Z2Ai applies on 

(E,J„,(A,),A), (44) 

where Ai = f[Ai]{Jni{0)), it will make the canonical form transform into 

{E,JnJ—^),A':), (45) 
\1 + Z2aioJ 

where A" = fl^^EA)A] ° f[(EA)j] ° "^"i ( i+zga o ) ' "^^^ entanglement classes parameterized by 
a = (Aj,aj„) are not inequivalent with different parametric values, while the parameters 
undergo a transformation characterized by Z2, i.e., 

I 

^ ) 1 + Z2aio _ ^^g^ 



A, 



a" 

"in 



Here, a"^ are elements of A". 
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3.3 The third superposition T^J'J-^ 

In this case 

Jn. + ^3A,=^f;i}(J„,(0)). (47) 

From the above equation it can be inferred that there should exist a similarity transfor- 
mation 

S-'glTpn.mS = JuM + zsa,o) • (48) 



Then 



S-'Jr^MS = gljtVnA^ + ^3aio)] (49) 



and 



= flAi] o glTJ]~^ o J^. {\ + zsttio) . (50) 
We note that when the T transformation of + z^Ai applies on 

(£;,j„,(A,),A) , (51) 

where Ai — f[Ai]{Jni{0)), it will make the canonical form transform into 

{E,JnA>^i + Z3aio),A'l') , (52) 

where A'-' = /[^.j o^fj^,^] o J„. (A + zsttio). The entanglement classes parameterized by a = 
(Aj, ttin) are not completely incquivalent with different parametric values, while parameters 
undergo a transformation characterized by z^, i.e., 

' ^ . (53) 



/// 

"in 



Here, a"^ are elements of A'l 
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3.4 Results 



As far we know that different A are inequivalent up to the following transformations 



A, 

, 1 + Z2aio 



- f[{EA)A] ° f[{EA)J] ° ( 1 , ^ „ ) ' (55) 



A = f[A,]og'J'J-'oJ^^{\ + zsa,o) . (56) 
or their combinations in proper orders. 



(57) 



where V = llzT^z, A = ©A- 



Notice that the above-explained classification scheme should be applicable to the 2 x 
N X N case, a special case of the general L x N x N system. By then there will be only 
the superposition Tj^Jj], and {E, J) is just the entanglement class up to the symmetry [3] 



which is a special case of Eq. ([35 



3.5 Examples 

In the following we employ the new method to characterize the entanglement classes 
of2x N X N and 3 x N x N state under SLOCC. First, we show how the present method 
works by comparing it with the previous result [4] for a 2 x x system (Appendix A 
in [4]); then we use the present method to construct the canonical form for the specific 
3 X N X N system discussed above. 



3.5.1 2 X A^ X A^ case 



The quantum state of 2 x N x N can be represented by a matrix pair = (Fi, F2). 
The invertible operations (SLOCC operations) act on the three partite in the following 
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forms: 



(59) 



where T is a 2 x 2 matrix. Suppose one typical canonical form of an entanglement class 
is 

//••- \ /••. \\ 



vj/ = {E, J) 



\ 



1 
1 
1 



v 



■••/ 



A. 1 
A, 1 
A,- 



V 



(60) 



Because of the requirement r{E) > r{J), the admissible T matrix that keeps the rank 
inequality should be in the form of an upper triangular (^l^ ] [3J. For the sake of 
convenience we take the pair (ri,r2) to represent the two 3x3 submatrices hereafter. 
After the T transformation 



^11 + tuXi ti2 

r'l = I tu + ti2Xi ti2 I , r'a 

tll+tl2Ai 



the canonical form corresponding to it is now (i.e., P = Q T'^ 



^22^ t22 

^22 A j t22 

t22Xi 

l-p/-l\ 



(61) 



pr[Q,pr'^Q) = iE,j') 

( (■■. \ ( 

1 

1 

1 



•7 



V 



iii+ti2 





1 



tll+ti2Xi 







1 

tll+tl2Ai 
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(62) 



One can regard this as a symmetry among parameters Aj 



A' 



t22 Aj 



induced by 



tll+tl2Ai 

the transformation T in canonical forms. In the Appendix A of [4J, this symmetry is 
expressed as 



{E, J) ^ {E, J') 



(63) 



where the parameters in {E, J') are A- 



^22 Ai 



and thus the canonical forms with 



different values of the parameter actually belong to the same SLOCC class, that is, \i ~ A-. 
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By employing the method developed in this work, the quantum state can be manipu- 
lated as follows: 



^ = (ri,r2) = (i,A, + x) , (64) 

where x represents Jni(O). The transformation T makes 

{T[, r'a) = (til + tl2K + tl2X , t22K + t22x) . (65) 

The canonical form now can be obtained just by dividing the polynomial tu + ti2Ai + ti2X 
[the same role as T'f-^ before Eq. (l62l) ]: 



Here PT2Q can be expanded as 



^22^1 + t22X t22\ , ^11^22 ^11^12^22 2 fn^\ 

~ + ~ — rz — TT^^ ~ v. — ~ — m ^ ■ ' ) 



tu + ti2Ai + ti2a; til + ti2Ai {tu + ti2Ai) 2 {tu + ti2Ai 

Note that here J„;(0)^ = 0. The the Jordan form of the second matrix in Eq.f l66p is 
, , + Jr. (0), and the canonical form is then 

^^ = (1, A- + x)=(l, +x) . (68) 

V 111 + I^12\ J 

Compared to Eq. (!6^ we see that different parameters Aj and A^ = i^^^^^^x- belong to the 
same entanglement class induced by the transformation T, which is in agreement with 
the resuh of Eq.(l63]). 



3.5.2 3 X X iV case 



We consider the case of -i/^ = {A\^A2^A-^ where r{A\) = N. Here, the i/j can always 
be transformed into -0' = {E, A2, A'^) by invertible operators P,Q. The actual problem 
discussed in Sec. [3] is to construct the canonical form of {E, A2, A'^) under T, P, Q, where 
[^42, A3] = 0. With our method, this can be decomposed into two tasks: (1) the si- 
multaneous similarity transformation of matrix pairs {A2,A'^) (canonical form) and (2) 
application of superpositions induced by T among (£',^2,74.3) (symmetry between the 
canonical form). 
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In the triple- matrix form {E, A2, A^), its canonical form under P, Q is {E, J, A). The 
effective T transforms {E, J, A) into {E + ziJ + Z2A, J + z^A, A) whose canonical form is 
{E, J', A'). That is, T induces a symmetry (equivalent relation) between different canon- 
ical forms {E, J, A) and {E, J', A'). Our task is to find the relation between parameters 
in J', A' and J, A. The initial canonical form can be represented by a polynomial form, 
i.e., 



ij = {E,J,A) 
//•• 



v 



1 
1 
1 



v 



■••/ 



A. 1 
A, 1 
A, 



v 



■••/ 



CLii ai2 
ttjo CLii 

ttio 



•7 



= (1 , Aj + X , Ojo + ciiix + ai2X^) 
where x represents Jni(O). 



(69) 



In the following we show that the transformation T^^j^ in Sec. 13.11 may transform the 
above quantum state into the following form: 



(1 + ZiXi + ZiX , Aj + X , flio + CLil^ + (ii2x'^) 



(70) 



We transform the first matrix into a unit matrix by dividing the polynomial 1 + ziXi + zix 



and obtain the polynomial functions /j 

Xi + X 



) , 

[{EJ)jy 



x) and /, 



(21) 

[{EJ)A] 



[x). That is, 



aio + anx + ai2x' 

1 + ziXi + zix' 1 + ziXi + zix 



1, 



A, 



X 



ZiX 



l + z,Xi (l + 2iAi)2 (l + ziXi)^'' 
a^o (ctii — CLioZi + aiiZiXi)x 
1 + z,Xi (l + ^iA,)2 

0^12(1 + ^lAj)^ — 2l(«il — fljO^l + 0.ilZlXi) 



(l + ^iA.)3 



(71) 



= ' /Kij)J](^) ' f[{EJ)A]^ 

where only up to x"^ is needed due to J^.(O) = 0. And, then we transform the second 
matrix into Jordan form, which means the similarity transformation 



[{EJ)J\ 



'x) = J„ 



. A, 
T + ziXi 



(72) 
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In polynomial form it reads 



+ n , . ^ - M , . ^ ^, = ^nX , ,.x )= .,.x +^. (73) 



l + Z^Xi {l + Zi\Y (l + ^lA.)3 "'^1 + ^iA/ l + ^iAi 

with y = MxM~^, where x = implies y = 0. Hence, the quantum state is now in the 
form 

correspondingly. One can solve y from Eq.( l73|) by expressing y = Yl'ili like 

iEmiJnXYy^)) = ?/ = (! + ^lA.)' x + z,{l + z,X,) 'x' + -- - . (75) 
Take y into f\^E.j)A]^y) ^"^^ keep those terms up to x^, for the third matrix A[ we have 

+ (ttii - aio2;i + aii2;iAi) X + (1 + -ZiAi) x . (76) 



1 + Zi\ 

Finally, by comparing coefficients of different powers of x, we have 
(Aj, Ojo, Oil, ai2) ~ (Ai, a^O' '^ii' ^ 



1 + ziXi 1 + ziXi 



a.) 

'an-aioZi + anZiXi) ,ai2{l + ziXi)^ ] . (77) 



Similarly, we can get the result for T^^^J-^, i.e., 

(Aj , aio , an , 0^2) ~ (A- , a-Q , a-^ , 0-2) 

Aj flii Oj2 (1 + 0^20^2) ' 



1 + 220*0 ' 1 + Z2aio ' 1 + aioZ2 - aiiZ2Xi ' (1 + 22 (fljo - «jiAi))' 



3 \ ' 



(7^ 



and j 



(Aj , ttjo , «ji , ^12) ~ (A^ , a^o , flii , 0^2) 

l + anz^ (l + 23«*i) y 
More specifically, take Eq.( l79|l as an example. If a quantum state in the canonical form 
has one block with parameters (Aj , , an , 0^2) = (1,0,2,3), then it is SLOCC equiv- 
alent to the canonical form which has the corresponding block of (Aj , a'^Q , a-^ , 0-2) = 
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5' (T+fli)^)- order to complete the classification, in some special cases, i.e., 
when rescaling the whole matrix of J or A as discussed before Eg. (1221) . zero appears in 
denominators in the fractions of Eqs. fl78l) and (179|) . and the states have to be treated sep- 
arately. In all these cases, the method of similarity transformation plus the polynomial 
form of symmetries shall be applied iteratively. 

4 Conclusion 

In this work, we propose a general classification scheme for entangled quantum sys- 
tem L X N X N. In this scenario, based on the commutativity of the ILOs T, P, Q, 
the classification procedure is decomposed into two steps: the simultaneous similarity 
transformation of a commuting matrix pair into a canonical form and application of the 
internal symmetry of parameters in the canonical form. This is innovative to the gen- 
eral entanglement classifications, because by this scheme we can always extract symmetry 
properties from the general equivalent relation and then leave a relative simple canonical 
form for the entanglement class. In most cases the major challenge comes from the de- 
tailed forms of the representation of symmetries. For demonstration, a concrete example 
of entanglement classification for a type of 3 x x system is presented, of which the 
symmetry properties are expressed in forms of polynomial functions. 

Acknowledgments 

This work was supported in part by the National Natural Science Foundation of 
China(NSFC) and by the CAS Key Projects KJCX2-yw-N29 and H92A0200S2. We thank 
Bin Liu for initial collaboration on this work. 



19 



Appendix 



A The structure of matrix A 



I. If the Jordan form J = ©J„-(Ai) is nonderogatory, that is, every Aj has geometric 
multiphcity 1, then matrices commuting with it can be expressed as A = (BiAi, where Ai 
is in the form of Eg. ( ITT]) (see, for example, theorem S2.2 in [TO]). 



II. If some of the Jordan blocks of J have the same eigenvalue, 

J = ■ ■ ■ © Jni(Al) © J„2(A2) © ■ ■ ■ © J„,(Ai) © ■ 



^0) 



Ai=A2=---=Ai 

where rii > rij for i > j. (We can always gather the Jordan blocks with same Aj together.) 
The matrix commuting with this block can be expressed as 

/An Ai2 ■■■ Aii\ 

A21 A22 ■■■ A21 
A = . 

\Aii A12 ■■■ All J 

where A^ is an x rij upper triangular Toeplitz matrix. As the main subject of our 
work is the study of parametric symmetries, we refer to Eqs. (S2.13— S2.15) in [10] for 
the detailed definition of Aij. Here we give an example of rii = 3, 77,2 = 2: 

/ O-llO O-Ul 0-112 0,120 Oi2l^ 

aiio fliii ai2o 

aiio . (81) 

a210 O211 0220 CL221 

\ a2io 0220/ 
This can be represented phenomenologically by the following matrix of polynomials: 



All A12 
A21 A22 



/[Aii](a^) /[Ai2](-^) 



Olio + OiiiX + aii2X 



ai20 + ai2ix 

^220 + «22ia; 



2) 



a2ioX + a2iix^ 

where x = J„i(0) given that rii > 72,2, and the last rii — n2 rows and columns of Eq. (!82|) 
should be omitted to give the form of Eg. (181 p . In the discussion of the parametric sym- 
metry of the canonical form in Sees. 13.21 and 13. 3[ the supposition and inverse functions 
applied in Eq.f l82|) can be proceeded directly with the only difference being the emergence 
of the nondiagonal elements, i.e., /[A12] ^"^^ /[A21] Eq.( l82|) . 
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